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Abstract 

In a SUSY GUT seesaw scenario, the largeness of the atmospheric neutrino mixing can reflect 
itself into an enhanced flavor changing mixing of beauty and strange right-handed scalar quarks. If 
the CP violating phase in such down type squark RR insertion is the main source of CP asymmetry 
in Bd fKs and the gluino contributions to K^-K^ and B^-B^ mixing are negligible, there is a 
correlation between Re {e Bk) and B^ —> fKs CP asymmetry, in addition to that with the strange 
quark CEDM. The current data on Re [e jcK) = (16.7 ± 2.6) x 10“^ imply that S^k should be 
greater than ~ 0.5 (0.25) for /xtan/3 = 1 (5) TeV, assuming the RR dominance in 6 —> s transition 
and the minimal supergravity type boundary conditions for soft parameters. 
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I. INTRODUCTION 


The flavor changing neutral current (FCNC) processes induced hy b —>■ s transitions have 
played a major role in probing proposals of new physics at the electroweak scale. While 
the study oi b sy keeps being the highlight in such investigations, more recently the 
relevance of 6 —> s induced purely hadronic decays has been emphasized, in particular in 
relation with the issue of testing CP violation in B physics. Among such decays, the process 
Bd —>■ (pKs has aroused much interest at least for two reasons: i) the process can occur 
only at the loop level in the Standard Model (SM), hence making it particularly suitable 
to spot sizeable contributions coming from new physics; ii) in the SM the CP asymmetry 
in such decay arises only from the indirect CP violation of the B mixing, hence one can 
safely state that within the SM the golden mode Bd —>■ J/'ll)Ks and the decay Bd —> (t)Ks 
yield the same amount of CP asymmetry, namely S^k = if the new physics entails the 
presence of direct CP asymmetry in the decay amplitude of Bd (t>Ks, this can be revealed 
by a departure from the equality between the two mentioned CP asymmetries Biiii- 
Indeed, the first data on the CP asymmetry in Bd —>■ (pKs iii 2002 gave rise to a wave 
of interest on this potentially very interesting signature of new physics which still continues 

today 0,i,[I3EIlIIlliiEi,IiSlil. 

Indeed, even though the discrepancy of the data with respect to the SM predictions has 
been constantly decreasing in the two years elapsed from the first results in 2002, still the 
current world averages of S^j^K and S^k are Q 

S^K = (0.34 ± 0.21), S^K = (0.726 ± 0.037), 

namely S^pK is about 2 a lower than the SM prediction S'™ ~ S'™. 

Notice that even if one does not wish to speculate much on the single above discrepancy, 
it is interesting to note that the data concerning b —>• qqs processes with q = u,d,s exhibit 
an overall discrepancy in the values of the CP asymmetry with respect to the SM predictions 
(indeed, the CP asymmetries in 5 —>■ <pKs and Bd ^ rj Ks turn out to be smaller than that 
in Bd J/'ijjKs which is measured from the tree le vel p rocess b —>■ ccs). This discrepancy 
might be a signal for physics beyond the SM [ll, [l^ HH . 

Assuming that the current low value of S^k is a signal of new physics, we need a new CP 
violating phase in b ^ s transition. An attractive possibility for such new physics beyond 
the SM occurs in supersymmetric grand unified theories (SUSY GUT’s) scenarios with 
seesaw mechanism for neutrino masses and mixings. In such scenarios, the large atmospheric 
neutrino oscillation can be related with a large 6 —*• s transition through down type squark 
and gluino loop effects. This flavor changing effect is parameterized by a mixing parameter 
i^ 23 )RR with a CP phase ~ 0(1). We follow Ref. for the definitions of the mass insertion 
parameters ((5 )^)^b’s. Let us note that the following relations 

(Sf^RR = {S^Yrr, (Sf^LR = (5?.)k, 

will be useful, when we relate {S^Yrr and {b^ 2 )RR ^^is work. For low tan/3, the single RR 
insertion can lead to some deviation in Sp^K, if gluinos and squarks are relatively light. For 
large tan/3 case, the double mass insertion can lead to an effective RL insertion of 0(10“^), 
leading to a significant deviation in S^px from the SM prediction. In Fig. 1 (a), we show the 
Feynman diagram for b ^ sg involving a CP violating {S^Yrr- 

The new CP violating phase in the RR insertion can affect the strange quark chromo¬ 
electric dipole moment (CEDM) through triple mass insertions, if there is an LL insertion 
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between third and second generation down squarks. [Fig. 1 (b) ] B H13. Since 
the LL insertion is generically present in the minimal supergravity (mSUGRA) case, the 
strange quark CEDM puts a strong constraint on the possible deviation of S^k from the SM 
prediction. However, the substantial theoretical uncertainties occurring when one relates the 
quark CEDM’s with the hadronic EDM’s suggest that it would be preferable to have some 
other observable at disposal in addition to the strange quark CEDM in order to constrain 


S(I)K- 

In this letter, we point out that the phase in the ((523)r?J? mixing parameter that would 
affect S^K can also contribute to direct CP violation within the neutral kaon system, namely 
Re (e /ex) through triple mass insertion. The Feynman diagram for Re (e /ex) with triple 
mass insertion [ Fig. 1 (c) ] is very similar to the Feynman diagram for the strange quark 
CEDM [ Fig. 1 (b) ]. Needless to say, making use of the observable Re (e /ex) to constrain 
some SUSY soft breaking parameters also entails theoretical uncertainties mainly ascribed 
to our ignorance in the evaluation of the relevant hadronic matrix elements. Our discussion 
shows that, even taking into account such a huge degree of uncertainty. Re (e /ex) still 
constitutes a precious tool in constraining the interesting flavor changing mass insertion 
parameter (^ 23 )^^^ ^^^1, in any case, it plays at least a complementary role to the strange 
quark CEDM in performing such a task. Although the large RR mixing is well motivated in 
SUSY CUT plus seesaw mechanism, we should emphasize that the link between Re (e /ex) 
and S^x is a generic feature of RR dominance scenario which in general can arise in other 
contexts. 

For definiteness, we will work in the mSUCRA boundary condition at the reduced Planck 
scale M* ~ 2.4 x 10^® CeV, namely the flavor universal scalar mass parameter mo, and the 
trilinear couplings A which is proportional to the Yukawa couplings. Then flavor changing 
off diagonal squark masses will be induced by the renormalization group (RC) evolution 
from M* down to Mqut — 2 x 10^® CeV, and subsequently from there to the Mw scale 0: 


{m\)ij = -^[V^XlV]ij(?>ml + A^) [Slog 


Stt^ 

from which one can estimate IlSl If 


M* 


M 


+ log 


M, 


GUT 


GUT 


Mw ) ’ 


( 1 ) 


(<5i3)ll ^ 8 X 10 ^ X e 


,-i2.7 


The LR or RL mixing will be small in this limit, and we ignore them in the following. 

For the RR insertion, it is also generated by the RC evolution. In the presence of a 
seesaw mechanism to give rise to neutrino masses, there will be the new Yukawa couplings 
responsible for the Dirac entries in the neutrino mass matrix. In a unihed context, such 
couplings intervene also in the evolution of the right-handed squark masses in the interval 
from the scale of appearance of the soft breaking terms down to the scale of right-handed 
neutrino masses. If at least some of such couplings are large (for instance one of the Dirac 
neutrino couplings could be of the order of the top Yukawa couplings), then one can expect 
sizeable contributions to the RC induced off-diagonal entries in right-handed squark mass 
matrix. Let us work in the basis where the down quark and the charged lepton mass matrices 
are diagonal. In this particular basis, we get 




Stt^ 

1 

Stt^ 


[Y^Ytv]^( 3mo A^) log ■ 


M, 


GUT 




yljyL]ki[VL]kj{‘i'ml + A^) log 


M* 


M 


GUT 
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Here Yat denotes the 3x3 Yukawa matrix of the term inducing the neutrino Dirac mass, 
which is diagonalized in the form, 




Y]\i diag (jjuii yv 2 i l/t's)) 


Ql = diag 



(3) 


where is a unitary matrix with a single phase in the standard parameterization, [/at is a 
general unitary matrix, and the phases are subject to the constraint, + 02 ^^ + 03 ^^ = 0 . 
Hereafter, we assume that the right-handed neutrino mass matrix Mat is diagonal in the 
basis where Yjq is diagonal. 

Mm = uI^MmU*^, Mm = diag (Matj, Matj, • (4) 


The simplest case of this sort would arise for Mm proportional to a unit matrix. A case with 
non-degenerate eigenvalues is also possible. Such a texture with simultaneously diagonaliz- 
able Ym and Mm may result from simple U(l) family symmetries. Under this assumption, 
coincides with hermitian conjugate of the usual definition of PMNS lepton mixing matrix, 
which is not always the case for arbitrary Mm- 

From Eq. (( 21 ) , we can estimate In particular, having in mind a large Yukawa 

coupling of the third generation, the entry ( 5 f 3 )i?_R would be largely affected: 


(<^23)^2^ ^ 2 X 10 ^ 


( Mm, \ 
VIO14 GeVy ■ 


(5) 


Its size depends on the right-handed neutrino mass scale, because in Eq. (j2I) grows with 
Mm 3 for a fixed value of neutrino mass. Sizes of (^^ 2 )^?^? ('^ 13 )^^ crucially depend on y^,^, 

yu 2 , and [Ul] 3 i as well. In general, they can be large enough to influence ex, Re (e /ex), and 
mixing signihcantly. Suppose that the neutrino mass spectrum has normal hierarchy 
with negligible lightest neutrino mass. In this case, Eq. ((2)) implies that {5f2)RR and 
scale as 


('^12)rr| ^ l(y3)RR| max |[Vc]3l|) X |(i523)rr 


( 6 ) 


where we regard all the elements of Vl as 0(1) except for [Ul] 3 i. Maximal size of ( 523 )flii I^at 
can be expected from Eq. m is 0(0.1) for neutrino Yukawa couplings to be perturbative. If 
the right-handed neutrino masses are degenerate so that yl^/yl^ ~ 0 . 2 , the resulting {5f2)RR 

2 ^ . Also I [Ul] 3 i| close to the current upper bound 


can result in a huge contribution to ex [ 2 / 


~ 0.15 j26|] can give rise to a similar size. For large tan0, this can lead to contribution to 
Re {e /ex) as well, which is bigger than its experimental value |27l. l28l |. The other mass 
insertion may modify B^-B^ mixing to a sizeable extent. 

One can make an observation from Eq. Q that (^ 12 )^^ and {5f^)RR will be highly sup¬ 
pressed provided that neutrino Yukawa couplings have strong hierarchy and that [Ul] 3 i is 
vanishingly small. In fact, the former condition is naturally realized in a scenario with 
SO(IO) unihcation j^. In this scenario, the neutrino Yukawa couplings are unihed with the 
up-type quark Yukawa couplings, thereby resulting in yl^/yl^ ~ lO”"^. Since y‘l^ly‘1^ is much 
smaller than eigenvalues of Mm should be split accordingly for the seesaw formula 

to yield correct neutrino mass spectrum. Then the RG induced (^ 12 )/?/? and will 

be small, and SUSY contributions to Re {e jex), and and B^-B^ mixings can be 

safely neglected. We will consider the correlation between S^px and Re (e /ex) in such a 
case, r eleg ating the more general case with sizable { 5 / 2 )rr and {Sf^)RR insertions to a future 
study j29|. 
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This paper is organized as follows. In Sec. II, we give basic formulae for the SUSY con¬ 
tributions for Ai? = 1 effective Hamiltonian which are relevant in the mSUGRA boundary 
conditions, and discuss how to include hadronic uncertainties in S^j^K- In Sec. Ill, we give the 
relevant information on AS* = 1 effective Hamiltonian and hadronic uncertainties related 
with Re (e /ex)- The numerical analysis is given in Sec. IV, and the results are summarized 
in Sec. V. 


II. CP ASYMMETRY IN Bd 4>Ks 


Let us start with Bd <pKs decay by recapitulating the effective Hamiltonian for AB = 1 
relevant to Bd —>■ (pKs With the operator basis in Eqs. (9) of Ref. P, it is given by 


H. 


AB=l 

eS 


r r 

^ + C2OI + Y, Q(/^)0.(/i) + ^7,07, + 

^ p=u,c L ^=3 


+ H.C. (7) 


One also has tilded operators Oj= 3 ,,,,, 10 , 77 , 89 , which are obtained from OiS by making chirality 
flip L ^ R. Then the Wilson coefficients for the tilded (chromo)magnetic operators, up to 
the second order of mass insertion approximation, [ Fig. 1 (a) ], read 


X Gp 77 _ TTOs 

V2 

_ 

mb 


A A nrrt _ 

-gMii^){ 62 ARR + g^M2{x){62ARRi533)RL 

Um3(x) - ^M4x)) 

Am2(x) + (4)kh(« 


( 8 ) 


' 33 )^ 1 - 


The loop functions for single and double mass insertions can be found in Refs. j2^ and 
[ 2 ^ . respectively. We also include the contributions from U 3 ,..., 6 . We have ignored the terms 
depending on { 623 )rl in the tilded Wilson coefficients, and and ( 60 A t.r i n the untilded 

Wilson coefficients, because they are all small within mSUGRA scenarios 1241 . 

We calculate S(j,K using the QCD factorization in the BBNS approach j30|. There are 
theoretical uncertainties from the divergent integral in the hard-scattering (H) and the weak 
annihilation (A) contributions such as dy/y. We adopt the suggestion by BBNS as follows 
0: 

[ dy/y ^ (1 + \og{mB/\h) 

Jo 

with 

0 < Pa < 1, 0 < (p77, ifA < 27r. (9) 

This prescription is an intrinsic limitation of the BBNS approach, and the associated un¬ 
certainties cannot be reduced at the moment. It turns out that these uncertainties are not 
very large if squarks and gluinos are relatively heavy, 350GeV < m, m^, but can be large 
for lighter squarks and gluinos close to the current lower bounds. If we assume the gaugino 
mass unihcation within SUSY GUT’s, the LEP bound on chargino {m^+ > 94 GeV) implies 
that rrig > 400 GeV. In the numerical analysis, we use rrig = m = 500 GeV, and the hadronic 
uncertainties become smaller. 
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Finally, we will impose the usual bounds that any new physics involving & —> s transitions 
should satisfy 

2.0 X 10"^ < BiB Xs-i) < 4.5 x 10“^ 

r r (10) 

AM, > 14.9 ps"h 

We should comment on constraints from b —*• another important ingredient. The 

branching ratio of the exclusive decay mode B —*• Kt^t~ is 

^ iF£+r) = (5.4 ±0.8) X 10"^ (11) 

which is consistent with the SM prediction j^. A model independent analysis in Ref. [s^ 
has shown that the region on the plane of (C^^, allowed by this decay at the 90% 
C.L. is an annulus with radius ~ 7 and thickness ~ 5 including the origin, where and 
are new physics contributions to the Wilson coefficients of the 5 —> s l'^l~ four fermion 
operators under their convention. On the other hand, for an 0(1) value of \{ 323 )ll\ and 
rrig = rh = 500 GeV, a maximal gluino-squark loop contribution to has the size ~ 0.2, 
which is much smaller than the extent of the allowed annulus. Although the data quoted 
in the above reference is less precise than the present one, the overall feature should remain 
the same. Also, incorporation of the inclusive mode does not make much difference. From 
this we can deduce that the effect of an RR insertion, which we are considering in this work, 
should be equally insignihcant. As for the effective RL insertion, it is much more severely 
constrained by i? —> Xgj than B —> Xsi^i~. For these reasons, we do not include the 
b —>■ s constraints explicitly in our analysis. 


III. Re {e /ex) 


The AS* = 1 effective Hamiltonian is given by 
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ttAS=1 

-^eff 


— GiOi ± C 2 O 2 ± CiOi ± ± CggOgg ± h.c. 


( 12 ) 


i=3 


where the operators O-r-y and Ogg are the same as in Eq. (|7j) except the replacements 
(s, b) —(d, s). The leading contributions to the Wilson coefficients of the (chromo)magnetic 
operators are provided by the triple mass insertions [ Fig. 1 (c) ], reading as 


C7y = 

Cs,= 


ttq;, 4 m 




('^13)^1(^33)111(^32) nil, 


9 m 


TTQ!, TUg 

ni^ m. 


^Ai(a:) ± ^^2(0;) 


('^13)LL(<^33)Lfi(<^32)RR- 


(13) 


The loop functions Ni[x) and N 2 {x) are available in Ref. 113 In these expressions, we have 
omitted double and triple insertion terms from other mass insertion parameters as they are 
small compared to the terms we kept above. As in Cjy and Csg, we ignore the (di 2 )LL 
contributions in 03 ^,,,^. The importance of the single LR insertion and the double insertion 
contributions to Re (e /ex) was pointed out in Refs, and 00, respectively, in the 
general MSSM frameworks without relying on the flavor universal boundary conditions. In 
the following, we show that the triple mass insertion can give an important contribution to 
Re {e /ex)- 
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After accomplishing the computation of the relevant SUSY contributions in the effective 
hamiltonian responsible for Re (e we now turn to the delicate issue of the hadronic 

uncertainties in the evaluation of Re (e /ei^). For dehniteness we will follow the treatment 
provided in Refs, js^, |3^. In the SM contribution, the main uncertainties reside in the 
evaluation of the B parameters and in the estimate of the strange quark mass 

mg. We dehne the non-perturbative parameters Rq and Rg as 


- r(1/2) 


-^6 = Bq 


121 MeV 


1 2 


_ms{mc) + md{mc) 




121 MeV 


n 2 


_ms{mc) + md{mc 


(14) 


The value of is rather well estimated both from lattice QCD and from analytic 

non-perturbative approaches [s^. In what follows, we employ the range of Rg, 


Rg = 1.0±0.2. (15) 

On the other hand, the situation of is very unclear, and there exist results from 

different approaches ranging within a factor of 2.2. For instance, the large-W limit predicts 
Re = 1. It is difficult to attribute a reliable uncertainty to such estimates. Taking Re = 
Rg = 1.0 results in the SM prediction of Re (e /e^^), which is smaller than the observed value 
by about 3 x 10“"^. According to Ref. j^, the best ht from Re {e /ex) yields Re = 1.23. 
However, let us emphasize that this is the case within the SM, and we cannot rely on a SM 
ht here in the presence of new physics affecting kaon decays. As a representative value, we 
use 

Re = 1.0 ±0.2. (16) 

In any event, our conclusion does not depend signihcantly on the precise value of Re or Rg 
because Re (e /ex) is dominated by chromomagnetic contribution in most of the parameter 
space. 

As for the strange quark mass, we make use of the range: 


ms{mc) = 115 ± 20 MeV. 


(17) 


For completeness, we also specify our choice of the isospin breaking parameter Hib. We use 
the value 

fliB = 0.06. (18) 

The uncertainty of Hjb ~ 0.08 could shift Re {e jex) by about 0.8 x 10“^, and neglecting it 
does not affect our conclusion. 

It is time to come to the main uncertainty in the SUSY contribution to Re (e /ei^). This 
is related to the evaluation of the matrix element: 


{Qa )o - A/ o 


3 11 {qq) 


2167r2 R3 


ml Be, 


where 




Am. 


Ogg — O, 


85 


(19) 


( 20 ) 


The uncertainty in the above evaluation is encoded in the value of the parameter Bq. The 
result of Ref. corresponds to Bq = 1. Unfortunately lattice computations are still unable 
to come up with a reliable estimate of the matrix element {Q~)o] in fact, even the sign of 
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this parameter is not certain yet, although the above reference estimated it to be positive. 
If we assume the opposite sign of Bq with the same magnitude, the SUSY contribution flips 
its sign. In spite of all this and even allowing for an uncertainty of a factor 4 in the estimate 
of \Bg\ = 1-4, we will show that the constraint on {S 23 )rr from Re (e /ex) still remains 
meaningful. 


IV. NUMERICAL ANALYSES 

Now it is straightforward to calculate Re {e /ex) and S(j,x from s ^ dg and b —>• sg, when 
the (^ 23 )^^ is the main new physics contribution beyond the SM contributions. Here we 
assume that SUSY contributions from {Sf 2 )RR and { 6 f^)RR are negligible as is the case for 
hierarchical neutrino Yukawa couplings and vanishing [Rijai. If we relax this assumption, we 
should regard 7 as a free variable that can be varied within a range compatible with other 
data such as B —>■ X^'j- Then, the calculated value of Re (e /ex) will change by a fraction 
of its SM prediction. Nevertheless, the qualitative feature remains true that the size of Re 
(e /ex) bounds the deviation of S^x from S^. 

We would like to point out the main result of this work in a simple way, before we give a 
detailed numerical analysis. If the {d'/^^)RR mixing is the dominant new physics contribution 
to Bd —>■ (j)Ks, we find the following from Eqs. m and (jHl) in the previous sections : 

Re(e/e^) : Cf7''(AS = 1) a /,(x) (21) 

S4,K : C|7''(AB = 1) OC f2(x) (73)hh + Mx) {il3)RL{ii)RR, (22) 

nib 

where /i=i, 2 , 3 (x) are the loop functions obtained in the previous sections. Now, if the SUSY 
contribution saturates Re (e /ex), then it is well known that one has to satisfy 

\{Si3)ll{5^3)lr{5^2)rr\ ^ 10 ”^ 

with an 0(1) phase ji^. Since the RG evolution generates {Sf^)LL ~ within mSUGRA 
scenario, we can derive the following upper bound: 

l('533)i^?('^32)Rij| ~ 10”^- (23) 

Note that this combination enters the calculation of S^x and B —through 
Csg[ 7 ^){AB = 1) along with { 6 ^ 2 )RR- For a small /itan/3 (corresponding to a small (^fa)/?!,), 
one can have larger {d^ 2 )RRj which is constrained by the lower bound on AMg and the 
B —>• Xgj branching ratio. For a large /itan/3 (corresponding to a large (d 23 )RL), (^ 32 )/?^ 
should be smaller in order to satisfy (El. In either case, we can expect that the deviation 
in cannot be that large for such {S^ 2 )rr satisfying the Re (e /ex) constraint, (El- 
Having described the qualitative features of our main points, we now provide the detailed 
analysis including theoretical uncertainties in S^x and Re (e /ex) as summarized in Secs. H 
and HI. We use the parameterization, (( 523 )iiJ? = We fix the modulus r at a maximal 

value compatible with B{B ^ ^sl), and vary its phase 0 from 0 to 2tt. For ptan/3 = 1 TeV, 
we use r = 1 which is dictated by requiring the validity of the mass insertion approximation. 
For /itan/3 = 5 TeV, we set r = 0.33, which is defined by the up^r bound on B{B —> X^q). 
Some values of <j) result in AMg smaller than the lower bound jof, and they are discarded. 
For each value of cj) consistent with the AMg constraint, we plot a point in the plane Re 


I^k) - S(f,K, following the procedures in Refs, and P,[lH, respectively. In Figs. 2 (a) 
and (b), we show the plots for ptan/3 = 1 and 5 TeV, respectively, with m = nig = 500 GeV. 
The thick vertical error bar shows the current data on S^k, and the two dashed vertical lines 
delimit the experimental value of Re {e /ex) j^ . 

Re (e'/ex) = (16.7 ± 2.6) x 10"^ (24) 

The full black box shows our estimates of and Re [e /ex) within the SM. Its width and 
height are the uncertainties in Re (e /ex) and S^x-, respectively. In Fig. 2 (a), we show a 
curve for /itan/? = 1 TeV, and r = 1 which may be regarded to define the boundary of a 
region of the two observables generically predicted in this scenario. For this curve, we £x 
Be = 2 and use the central values of Rq, Rg, and nig, as given in Sec. III. As mentioned 
previously, disconnected parts of the curve are excluded by the AMg constraint. If we 
turn on hadronic uncertainties, this curve gets broadened into the gray region around it. 
We estimate the uncertainty of an observable by taking its maximum and minimum values 
reached while varying the relevant input parameters in the ranges quoted in the previous two 
sections. For Re (e /ex), they are Rq, Rs, nig, and Be, and for they are qh, ^h, Qa, 
and ifA- Each of the uncertainties in Re (e /ex) and S^j^x is displayed by the horizontal or 
vertical error bars at hve selected points. The gray region is drawn by varying all the eight 
parameters simultaneously. It therefore covers wider space than is obtained by quadrature 
addition. Suppose that the sign of Bq is negative. The resulting curve and the region around 
it can easily be guessed by taking the mirror image of the present one around the vertical 
axis passing through the SM point. Even then the Re (e /ex) data gives a strong constraint 
on the possible value of Sf^x- 

We repeat the same exercise for /itan/3 = 5 TeV in Fig. 2 (b). We find that Sfpx > 0.25 
if the data on Re (e /ex) is imposed. 

Note that the constraint from Re (e /ex) is comparable to that from the strange quark 
CEDM. In particular the positive (negative) S^x is correlated with the positive (negative) 
Re (e /ex) within minimal SUGRA boundary conditions. Therefore, the old Belle data 
with the negative S^x implies a negative Re (e /ex) in the RR dominance scenario such 
as SUSY GUT models with right-handed neutrinos, which is clearly excluded by the data 
Re [e /ex) = (16.7 ± 2.6) x 10“^. If the old Belle data were still valid, then the RR 
dominance scenario should have been discarded. Our results provide a meaningful correlation 
between S^x and Re (e /ex) despite large hadronic uncertainties in both quantities. This 
is independent of the strange quark GEDM constraint, and probably has less theoretical 
uncertainties. 

One may wonder why we are not considering triple mass insertion on a squark line in the 
box diagram for ex- The size of the effective insertion, 

{^{2)fK ^ {^{^)lMz)lM2)rr. (25) 

is always smaller than 2x 10“"^ due to B{B —> Y^y). If we require that the SUSY contribution 
to Re (e /ex) be smaller than its experimental value, we get 

< 10 -', ( 26 ) 

which implies that 
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This limits the SUSY contribution to ex below 1/30 of its experimental value. In view of 
theoretical uncertainties in predicting ex, we may regard ex being always safe provided that 
Re (e /ex) constraint is satisfied. 

Let us add a remark on the CP violation in Bg-Bg mixing, whose indirect CP asymmetry 
nearly vanishes in the SM. The large RR mixing leads to a considerable modification in the 
mixing amplitude. This is evident from the fact that part of the curves in Figs. 2 was 
excluded by the AM^ constraint. Since this new piece of amplitude has a phase different 
from the SM one in general, the indirect CP asymmetry in Bg-Bg mixing can have a value of 
0(1) according to the phase of ((54)_Rii;. This will show up in the decay channel Bg —*• J/^d) 
as the time dependent CP asymmetry therein. For quantitative analyses, see Refs, 
for instance. 

In SU(5) SUSY GUT, the left-handed sleptons and the right-handed down-type squarks 
are tied to a 5, and therefore their flavor changing effects are related to each other. Off- 
diagonal elements of the left-handed slepton mass matrix are given by 

(m?),j ~ -^ylJ[VL]k^[Vl]kj{d>ml + A^) log (28) 

in a way similar to Eq. 0. In order to estimate lepton flavor violation from this mass 
matrix, we should go to the super CKM basis. If the down-type quark and the charged 
lepton Yukawa matrices are the same at the GUT scale as 

U = rt. (29) 


the right-handed down-type squark and the charged slepton mass insertion parameters are 
unified as well. Under this assumption, the current upper bound on the r —> py branching 
ratio ji^ . 

B{t /iy) < 6.8 X 10"®, (30) 


translates into the limit [1 




(31) 


for tan/3 = 10. This is roughly 1/10 of the size of {S 23 )rr we used for ptan/3 = 5 TeV, 
and we cannot expect considerable change in S'^x satisfying the r —py constraint. The 
assumption of Yukawa unification, however, leads to an incorrect mass relation. 


md rUe 
rUg 


(32) 


and Eq. (Ei should be modified to account for the mass ratios of the first and second 
generation fermions. Even in this case, a mass insertion involving a third generation is not 
much affected, and r —*■ /ly remains a strong constraint on S^k- 

In this work, we are considering SUSY GUT with right-handed neutrinos as an example 
of a scenario that gives rise to large RR mixing. This is why r —*■ /ly constrains S^k- Here, 
we would like to stress again that the interconnection between Re(e /ex) and S'^x, unlike 
r —> py, is a common consequence of large RR mixing, which is not specific to SUSY GUT. 
Suppose that there is a large RR mixing but no LL mixing in the squark sector at the 
reduced Planck scale due to a flavor symmetry and that we do not have a unified gauge 
group. Even in this case, one has a strong correlation between Re(e /ex) and while 
r —> py is unrelated to S^k- 
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If we considered more general scalar masses at M* with some flavor structures, then our 
results will be changed accordingly. The Wilson coefficients for Cgg’s for both AB{S) = 1 
have to include other mass insertion parameters such as ('^ 12 )^^; i^ 23 )LL, etc., which 

were neglected in Secs. II and III because they are small within the mSUGRA scenarios. 
Still we should make it sure that the new flavor physics that affects Si^k does not contribute 
to Re(e /ex) too much, and this could make a strong constraint on new sources of flavor 
and CP violation despite theoretical uncertainties in Re(e /ex)- 


V. CONCLUSIONS 

In summary, we showed that if the RR b ^ s transition is large with 0(1) phase, 
it can affect not only S^f^x through double mass insertion and the strange quark CEDM 
through triple mass insertion, but it also affects Re (e /ex)- The correlation between the 
two observables are strong despite large hadronic uncertainties in both observables, within 
mSUGRA boundary conditions with flavor universal scalar masses at M*. The current data 
on Re (e /ex) indicates that S^x should be in the range of 0.25-1.0, which is now in accord 
with the present world average of S^px- 
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FIG. 1: Feynman diagrams for (a) b ^ sg with double mass insertions, (b) strange quark CEDM 
with triple mass insertions, and (c) s ^ dg with triple mass insertions, involving {S 2 ^)rr as the 
dominant source of new CP violating parameter contributing to and Re {e lex). 
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(a) ^tan/3 = 1 TeV (b) ^tan/3 = 5 TeV 

FIG. 2: SfpK vs. Re {e lex) for (a) fitanP = 1 TeV and (b) fj,taiif5 = 5 TeV, with m = rrig = 500 
GeV. Experimental bounds on Re (e /ex) and S^px are depicted by the vertical dashed lines and 
the thick vertical error bar, respectively. Their SM predictions are marked by the black box, whose 
extent indicates their uncertainties. The black curve does not include hadronic uncertainties, and 
the gray region includes them. The respective uncertainties in Re {e jex) and S^px are shown by 
the horizontal and vertical error bars at some selected points. 
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